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SUMMARY 


This paper describes a computer program developed at the Langley Research 
Center to fit smooth surfaces to the component parts of three-dimensional air- 
craft configurations . The resulting equation definition of an aircraft 
numerical model is useful in obtaining continuous two-dimensional cross section 
plots in arbitrarily defined plajies, local tangents, enriched surface plots 
and other pertinent geometric information. 

At Langley Research Center and other government and industrial facilities, 
the geometry organization used as input to the program has become known as the 
Harris Wave Drag Geometry. An early version of this point siirface definition 
was introduced in a NASA TM X-9^7* (See ref. 1.) 


INTRODUCTION 


An important part of the design, analysis, and construction of aerospace 
vehicles, as well as ships, automobiles, etc., is the construction of surface 
models for the required study. Mathematical models are the most useful because 
of their versatility and exactness of definition. Many more different mathe- 
matical models may be analyzed by a computer program than would be possible by 
any other method within a given time frame. 


A set of planes defined by points to approximate a curved surface is the 
simplest type of mathematical model. A large number of points on the surface 
must be given to achieve an accurate definition using a set of planes. 


The surface fitting technique used in this paper is based on approximating 
complex surfaces by piecing together smaller surface "patches." The surface 
equations were conceived eind developed by Steven Coons. The term "Coons' 

Patch" is associated with this piecewise definition. (See ref. 2.) The ccam- 
puter program under discussion uses the bicubic form of the patch eq\xation. 

(See ref. 3-) This form of the patch equation requires 48 parameters consist-” 
ing of the coordinate position of the corners of the patch (12 parameters), 
the derivatives of the variables at the corners with respect to the parametric 
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lndei>endent variables (24 parameters), and the cross derivatives with respect 
to the parametric variables (12 parameters). 

Although patch surface definition has been applied in many ccanputer pro- 
grams, the process of determining the derivatives has been varied. One process 
has been to aijply phantom points near the corner points to approximate the 
derivatives by finite differencing. Many programs, including the one described 
herein, set the cross derivatives equal to zero. The LRC siirface program iises 
parametric cubic spline functions to approximate the derivatives of the space 
variables with respect to the parametric independent variables. 

The organization of the geometry input data to the program is identical 
to the geometry used for several standardized aerodynamic analysis con5)uter 
programs. (See ref. 4.) By using the same input format throughout the design 
and analysis procedure, the designer can be assured of a consistent numerical 
model. 

It should be pointed out that the program is not limited to aircraft con- 
figurations since a component part such as a fuseleige can be some other object. 
Also, it is a simple matter to substitute an alternate data input format if 
this might prove more convenient. 

Ciirrent program options that are available are to present on an XY 
plotter or on a cathode ray tube orthographic projections and cross section or 
contour plots. The patch equations are used to compute the required vectors, 

A complete description of the program and its usage is docianented in 
NASA TM X-3206. (See ref. 5.) 


SYMBOIS 


A,B,C,D 

a,b,c,d 

F 

G 

L 

M 

P 

S 

t 


parametric cubic spline coefficients 
plane equation coefficients 
boiandary matrix 

matrix equation of the solution of a patch and a plane 
chord length 

blending function matrix 

a vector whose components are functions of t 
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component of surface patch equation 
independent variable in cubic equation 


u,w 


independent variables in patch equation 
a vector whose components are functions of u and w 


V 


coordinates of a point 
pitch angle 
roll angle 
yaw angle 


PROBLEM DESCRIPTION AND METHOD OF SOLUTION 


A functional definition of an aircraft numerical model is useful for 
obtaining pertinent geometric information for use in aeronautical analysis 
programs or for building models for wind-tunnel testing. A computer program 
has been developed to provide this functional definition in surface patch 
equations . 

The organization of the input data for the program is broken down into 
various component parts of an aircraft configuration such as wing, fuselage, 
nacelle, vertical and horizontal tails . A typical input configuration is given 
in figure 1. Sufficient data points with no abrupt changes in curvature must 
be supplied as the modeling technique is designed to approximate a smooth 
surface . 

A surface patch is a portion of surface bounded by space curves and is a 
fionction of two variables, u and w, which can have only the values between 
0.0 and 1.0. It is necessary to supply parameters for each patch equation. 

A U X 4 boundary matrix is constructed for each component of the patch equa- 
tion (X,Y,Z). It is called a bo\indary matrix because the parameters are func- 
tions of the patch boundaries: corner points, corner derivatives and cross 

derivatives. The cross derivatives are set equal to zero in this application. 

The boundary matrix is then pre- and postmultiplied by constant -valued 
blending functions and stored for further use in this form. The blending 
functions provide a blending effect from one surface patch to an adjoining 
patch. Each surface patch is represented in exactly the same matrix form which 
proves to be very convenient for computer processing. The basic patch equations 
are given in fig\ire 2. 

The comer derivatives for the surface patch equations are obtained by 
tising a three-dimensional parametric cubic spline curve fit technique. The 
coordinate points describing the curve are expressed as cubic functions of one 
variable. Consecutive points in the u-direction are fitted in this manner 
and also in the w-direction. 

In the curve fit method, a series of adjacent polynomial segments between 


x,y,z 

e 

0 
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each given point is used to represent the curve. The adjacent cubic segments 
are related by setting the second derivatives equal at the common points. 

A matrix equation is developed where the slopes are the \inknowns. The 
second derivatives at the first and last points on a curve are made equal to 
zero since there are two less equations than unknowns. Since the matrix 
equation is tridiagonal, the Thomas Algorithm, which is equivalent to Gaussian 
Elimination without pivoting, is used for the solution. The Thomas Algorithm 
is a rapid, recursive method for solving a tridiagonal system of equations. 
Figure 5 shows the basic equations used. 

Parametric cubic splines were chosen because of their several advantages. 
Parametric curves are not sensitive to infinite slopes, spline curves approxi- 
mate smoothness satisfactorily, and cubics are the lowest order polynondal able 
to twist through space. 

To create orthographic projections of the given body surface, angles of 
yaw, pitch, and roll miist be supplied. The three-dimensional rotation 
equations are applied directly to the patch equations, where the rotated patch 
equation has only two components for the desired paper plane. 

The surface plots may be enriched to any degree. Two equations in w 
may be obtained by giving u a value between 0.0 8ind 1 . 0 . These eqxiations may 
be solved by giving w values between 0.0 and 1.0 for additional body surface 
lines in the w-direction. Similarly, two equations in u may be obtained by 
giving w a value between 0.0 and 1.0, and solving by varying u from 0.0 to 
1.0 for additional curves in the u-direction. 

A partial "hidden line" test may be used to delete points that face away 
from the viewer. A surface normal vector is computed from the rotated and 
projected patch equation at each requested point. The computed point is 
visible and plotted if the siurface vector is positive. If the surface vector 
is negative, the point is not plotted. Figure 4 gives the basic equations for 
the surface plots and figure 5 presents some examples. 

The simultaneous solution of the patch equations and the equation of a 
plane is used to provide cross section or contour plots. The cutting plane 
is defined by three points and may be at any orientation. The result is a 
4x4 matrix composed of constant elements and may be evaluated for an equation 
in the two variables u and w. The problem then becomes that of finding the 
zeros of a cubic polynomial in one of the parametric variables. Since there 
are three zeros, the program has logic in it to pick one of the points. If 
an error is made, it normally stands out immediately. The points defining the 
cxirve of intersection are then rotated and translated so that the intersection 
plane coincides with the YZ-plane of the paper. Figure 6 gives the basic 
equations for obtaining the cross section or contour plots and seme typical 
output plots are given in figure 7« 
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CONCLUDING REMARKS 


A coniputer program has been written that fits smooth surfaces to the com- 
ponent parts of an aircraft configuration. The resulting surface patch 
equations C8ui he useful anywhere an equation definition may be required. Pro- 
gram options currently provide plotting of the enriched surface at any desired 
viewing angle and the solution of a number of patches and a plane for generating 
cross section or contour plots through the given sxirfaces. 

Although the program was written to accept a point surface definition of 
an aircraft as input, it is a simple matter to substitute another input format 
to describe another object. 

The program has been used for on-line display on a cathode ray tube and 
to drive Gerber, Vaxian, and Calconap plotting devices. 
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Figure 2 . - Surface patch equation . 
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Y(t) = P(t) * At + Br + Ct + D 
Z(t) = 

Cubic Coefficients in Terms of End Points and Slopes (for one segment): 
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System of Equations to Solve for Unknown Slopes: 
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Figure 5- - Parametric cubic spline curve fitting. 
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Two Component Rotated Patch Equation: 


s; - ^ ( ie.ip ) + Sy ) 

Sy' " ^4 (f e.tf- » + Sy fj ( f e,(ir ) + ( ^,0,,^ ) 



Figure 4, - Rotated eind enriched patch. 
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Figure 5. - Continued. 
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Patch Equation: ^(u.w) U W* 

V(u,w)'USW* = Y(u,w) - US W* 

Z(u,w) U S W* 

Z 

Equation of a Plane : 

a x + by + cz-d = 0. 
then: 

U [a + b S + c S^] W* - d - 0. 


with: 

G-aS +bS +cS, 

X y z 

U G » d * [ u l] 
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Figure 6. - Cross sections or contours through a patch. 





Figure ?. - Continued. 
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